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We present a method for reducing the order of ordinary differential equations satisfying a given 
^v^j ~ scaling relation (Majorana scale-invariant equations). We also develop a variant of this method, 

■ aimed to reduce the degree of non-linearity of the lower-order equation. Some applications of these 

f"^ ' methods are carried out and, in particular, we show that second-order Emden-Fowler equations can 

£SJ , be transformed into first-order Abel equations. 

The work presented here is a generalization of a method used by Majorana in order to solve the 
O |' Thomas-Fermi equation. 
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I. INTRODUCTION 

> ■ 

In a recent paper m we have described a method, originally due to Majorana M, able to give the series solution of the 
. Thomas-Fermi equation (with appropriate boundary conditions) through only one quadrature. Such a method, giving 
a (semi-analytic) parametric solution of the considered equation, is based on a particular double change of variables 
which transforms the second order Thomas-Fermi equation into a first order equation, whose solution is then obtained 
£^ ■ by series expansion. 

Here we show that the transformation method, used by Majorana in that particular case, applies to a large class 
• of ordinary differential equations as well, and prove a simple but general theorem for reducing the order of these 
equations. 

The Majorana idea is a straigthforward generalization of known concepts and to show this we briefly recall, in the 
following section, some definitions and peculiarities of particular differential equations. In section III we then introduce 
a new class of differential equations and give the method for reducing the order of such equations. In section IV a 
variant of this method is presented and in section V some applications are reported which are particularly relevant in 
mathematical physics. 
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H | II. PRELIMINARIES 

Let us consider a general differential equation of order n in the independent variable x and dependent one y: 

F(x,y,y',y",...,yW) = , (1) 
where a prime ' denotes differentiation with respect to x. 

Eq. (|l|) is said to be an autonomous equation if the variable x does not appear explicitly: 

F (x, y, y', y", . . . , yW) = F (y, y \ . . . , . (2) 

In such a case, by changing the set of variables from (x, y(x)) to a novel one (?/, u(y)) through: 

y' = u(y) 

y = u (y)—^ — (3) 

ay 
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where u(y) is a given function of y, the considered differential equation can always be reduced to an equation of order 
n — 1 in the independent variable y and dependent one u 13] : 



G \ y ^Ty^--dy^) = ° ■ (4) 

A differential equation (Q) is, instead, equidimensional-in-x if it is invariant under the transformation x — > ax for any 
a ^ 0: 

F (ax, y, a^ 1 ?/', a" 2 y", . . . , a~V n) ) = F (x, y, y', y" , . . . , y^) . (5) 

This equation can be transformed || into an autonomous equation in the variables (z, y(z)): 

( dy d 2 y d n y\ 



by changing the independent variable: 



dz ' dz 2 ' ' dz™ 

x = e z 

d _ d 

dx dz 



Thus, cquidimcnsional-in-x equations of order n can always be reduced to differential equations of order n — 1. 
Scale-invariant equations satisfy the property: 



F 



(ax, a c y, a^y' \ a c ~ 2 y", . . . , a^V^) = F (x, y, y', y", . . . , . (8) 



(i.e. they are invariant for x — > ax, y — > a c y) for any value of a ^ and some value c, and they can be transformed 
into equidimensional-in-x equations ||: 

du d 2 u d n u 
dx ' dx 2 ' ' dx 



Glx^ix),— -^,...^) = , (9) 



with G obeying Eq. (|J), by performing the following change of the dependent variable: 

y(x) = x c u(x) . (10) 

Even in this case, scale- invariant equations of order ncan be thus reduced to equations of order n — 1. 

Finally, differential equations which are invariant under the transformation y — > ay for any a / are said to be 

equidimensional-in-y or homogeneous equations: 

F (x, ay, ay', ay", oyW) = F (x, y, y' , y" , y^) . (11) 

By changing the dependent variable through: 

y{x) = e»W , (12) 
it can be transformed into an equation of order n — 1 in the variables (x,u{x)) pJ: 

„ ( , . du d 2 u d n ~ 1 u\ „ , „. 

Thus, we know that the order of a given differential equation can always be reduced of one unit if this belongs to one 
of the four different classes mentioned above. 
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III. MAJORANA TRANSFORMATION 



For future convenience we now consider scale-invariant equations from a different point of view and introduce another 
class of differential equations. Equation ([!]) is said to be Majorana scale-invariant if it is invariant for x — * a c x, 
y — > ay for any a ^ and some value c: 

F (a c x, ay, a^y', a x ~ 2c y\ a^y^) = F (x, y, y', y", y^) . (14) 

It is easy to prove the following proposition: a Majorana scale-invariant equation of order n can always be reduced to 
a differential equation of order n — 1. In fact, by changing the role of the dependent and the independent variable 
x — > y, y — ► x Eq. (Q) can be transformed into: 

g(v,*<v)£,%,.--,%) =0 • (15) 



dy ' dy 2 ' ' dy n 



where now G satisfies the condition (g): 



. r _ 1 dx c _ 2 d 2 x d n x\ n( , , alx d 2 x d n x\ 

G(ay,a c x,a c L — ,a c 2 — ,...,a c n —\ = G[y,x(y), ) , (16 
ay dy z ay 11 ) \ ay dy £ dy n 



that is Eq. ( |15| ) is scale-invariant and the order can be reduced of one unit. 

The concept introduced above is not a really new one, since it is related to that of scale-invariant equations. However, 
the reformulation of the problem in these terms is useful for developing the method for the implementation of order 
reduction, which is a generalization of that used by Majorana in the framework of the Thomas- Fermi equation [jj. 
We now describe in detail such a method, which carries out the solution of Eq. ([!]), with F satisfying Eq. (Ei), as 
given in parametric form: 

x = x{t) 

(17) 

y = y(t) 

Let us assume that x in Eq. (|17]) depends on the parameter t through the function y(t) and, eventually, on t itself: 

x = x{t,y) = x(t,y{t)) = x(t) . (18) 

Since (x(t), y{t)) in Eq. (|i"7f) is a solution of the considered differential equation ([!]), supposed to be Majorana scale- 
invariant, Eqs. jl7|) and (|l8|) must satisfy the relation (|l4]), meaning that for any a / and a given value c we 
have: 

x(t,ay) = a c x(t,y) . (19) 
This implies that x(t, y) should be an homogeneous function of y: 

x(t,y) = x(t,l)y c = zy c , (20) 

where z = z(t) can be considered as an arbitrary but given function of the parameter t. Note that, in such a way, the 
only unknown function to be determined in order to satisfy Eq. (Q) is y(t), and the parametric solution (|l7j ) can be 
rewritten, after Eq. (^p|), as: 

= mv c {t) 

(21) 

= y(t) 

with z(t) an arbitrary but given function of t and c is determined by Eq. jl^). 
We have now to translate the differential equation ([!]) for y(x) into an equation for y(t) in (|lj). In the following, 
differentiation with respect to t will be denoted with a dot ', while a prime ' refers to differentiation with respect to x 
as above. The i-derivatives of x are, from Eq. (pl|), as follows: 
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z + cz- ) y c = Xl (t,y,y) 

v; 



z + 2cz — h cz 
V 



(c-1) 



X2 (t,y,y,y) 



(22) 



x n (t,y,y,y,. . . ,y} 



Using these expressions we can obtain the ^-derivatives of y, which are present in Eq. ([!]), in terms of t, y and its 
^derivatives: 

y' = - = yi(t,y,y) 

X 



., xy — xy , . ... 

y = — = y?{t,y,y,y) 



(23) 



V {n) = ■■■ = y n (t,y,y,y\..., y ) . 
The differential equation for y(t) is then obtained by substituting Eqs. (pl|), ( p2] ) and ( |2^ ) into Eq. (Q): 

P(t,y,y,y,...,y} = , (24) 

where: 

F (t,y,y,y,...,y^ = F (x{t,y),y,yi{t,y,y), . . . ,y n {t,y,y,y, ■ ■ ■ • (25) 

Note that, from Eq. (|2l]), the function x{t,y) is homogeneous (with respect to y) of degree c, while the functions 
yk(t, y, y, y\ ■ ■ ■ , V) are homogeneous of degree 1 — kc: 

k\ l-kr ' k 



y k \t,ay,ay,ay,...,ay j =a ' c y k \t,y,y,y, . . . , y j , (26) 

as required for the Majorana property ( p"4| ) to be satisfied. In particular from this we also deduce that the differential 
equation in ( p4] ) is equidimensional-in-y, since: 

F (t, ay, ay, ay,..., a yj = 

(x(t, ay), ay, y x {t, ay, ay),..., y n {t, ay, ay, ay, . . . , a yfj = 

(a c x{t,y),ay,a 1 - c y 1 (t,y,y),...,a 1 - nc y n {t,y,y,y\...,y)j = (27) 
{x{t,y),y,yi{t,y,y), . . . ,y n {t,y,y,y, . . .,y)J = 



F (t, y,y,y,..., tf) 



that is the function F satisfies the relation ([11]). We can then use the transformation in ([L2|) to reduce the order of 
the equation. More precisely we set: 



u(t)dt 

y(t) = eJ , (28) 



so that the i-derivatives of y(t) are as follows: 



y = uy = Ux{u)y 

jj = (u + u 2 ) y = u 2 (u,u) y 



. .. 

V = ... = u n [u,u,u,..., u y 



(29) 



The unknown function is now u{t) and the differential equation of order n — 1, obeyed by this quantity is obtained 
by substituting Eqs. (§|) into Eq. (||): 

F{t,y,u 1 y,u 2 y,---,u n y) = (30) 

or, by using the homogeneity of the function F (F(t, y, uiy, u 2 y, ■ ■ • , u n y) — F{t, 1, u\, u 2 , ■ ■ • , u n )), in all points where 
y(t) is different from zero we have: 

F(t,l,ui,u 2 ,...,u n ) = . (31) 
In terms of the initial function F in Eq. ([!]), by noting that: 

x(t, 1) = z(t) 

U\ u 

yi(t,l,Ul) = 7— r = — = V!{t,U) 

Xl(t, 1, Ml) Z + CUZ 



y 2 (t,l,U!,u 2 ) 



x 1 (t, 1, ui)u 2 — X 2 (t 7 1, til, u 2 )ui 

v 2 (t,u,u) 



iu — + (1 — 2c)iw 2 + c(l — c)zu 3 _ 



(32) 



(i + cuz) 3 



y n (t,l,ui,u 2 , . . . ,u n ) = ... = v n [t,u,u,u, u 



and using Eq. (p5|), we have the final equation for u(t): 

F (z(t), l,ui (i, it), v 2 (i, u,«), ... ,w„(t, u,tt, it, ( u ' )^ = . (33) 
Summarizing, the parametric solution of a Majorana scale- invariant differential equation of order n has the form as 



in Eq. fl2j\ ) with z(t) an arbitrary but given function of the parameter t and y{t) is written as in Eq. (28), where u(t) 

satisfies the differential equation fis3[ ) of order n — 1 and the functions Vk(t, u,u,u, . . . / u ^ ) are evaluated as in Eqs. 
HP- 
TV. THE METHOD OF THE AUXILIARY FUNCTION 

Depending also on the choice for the function z(t), it could happen that Eq. (|33|), although is of order n — 1, is too 
much hard to be solved. In some cases the following procedure can be used to reduce the degree of non-linearity of 
Eq. @. 

Let us perform a change of the dependent variable: 

u 

u — > v = . (34) 

z + cuz 

In this case we have: 
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v 2 = (1 - cvz) - + (1 - c) v 2 , (35) 



2 



and the function y(t) in Eq. (SWj is given now by: 



vz 

-dt 



y(t) = eJ I -cvz ) ( 36 ) 
where satisfies the differential equation of order n — 1: 

1, u, (l-cwz)^ + (l-c)i; 2 , ... ] =0 . (37) 



V. APPLICATIONS 

As an application of the Majorana method described above, let us consider the Emden-Fowler equation: 

y" = x a y b (38) 

(with a, b two real numbers), which is of particular interest in mathematical physics 0]. It satisfies the relation jl^ ) 
with: 

c = ^| , (39) 
a + 2 

so that the method may apply only for a ^ —2. In this case, Eq. ([33]) for the Emden-Fowler equation (§81) is: 

zit — zu + (1 — 2c)iu 2 + c(l — c)zu 3 a 



(z + cuz) 3 

with c given in Eq. (|39|). After some algebra we arrive at the following first-order equation for u(t) 

du 



(40) 



< 41 > 



where: 



= a(t) 


+ (3{t) u - 


- 7 (t) W 2 


a(t) = 


z a z 2 




m = 


3cz Q+1 z - 


z 

f - 




z 


7(t) = 


3c 2 z a + 2 - 


f 2c - 1 








6(t) = 


z 


' c(c-l) 



(42) 



With the Majorana method we have thus transformed the Emden-Fowler equation ( |3q ) into an Abel equation ( |4l| ) 
of the first kind. Depending on the particular problem to be solved, the final equation can be further simplified with 
an appropriate choice for z(t) which, however, cannot be chosen equal to a constant (in this case z ^ and Eq. ( |39| ) 
would not be a differential equation for u). 

A relevant case is that of Emden-Fowler equation with 6=1 for which, from Eq. ([39]), we have c = and Eq. (|4] 
reduces to a simpler Riccati equation: 



^ = z a z 2 +'iu-u 2 , (43) 
at z 



which for z(t) = t becomes: 



G 



^ = t a - u 2 . (44) 

Another interesting particular case is that of Thomas-Fermi equation, which is an Emden-Fowler equation with 
a = -1/2, b = 3/2: 

y3/2 

- w ' (45) 

The corresponding first-order equation ( [i"l| ) for u(t), choosing: 

z{t) = [12(1 -t)] 2/3 , (46) 

is: 

du 16 / 1 \ (1 \ 2 2 ,„ , o 

This equation was obtained by Majorana Q in studying the Thomas-Fermi equation. 

Emden-Fowler equations can also be analyzed by using the method of the auxiliary function outlined in the previous 
section. In this case, Eq. ( |37| ) for v(t) becomes: 

at 1 — cuz 

where c is given in Eq. (|39|) . Note that also in this case we cannot choose z(t) — constant (i ^ 0) from Eq. ([36]). 
Although Eqs. (48) and ( |4l| ) are different, in the particular case with 6 = 1, and thus c = 0, we again obtain a Riccati 
equation: 

J = z (z a - v 2 ) . (49) 

Instead, following the method of the auxiliary function with : 

z{t) = 12 2 / 3 i 2 , (50) 
the Thomas-Fermi equation can be transformed into the first-order equation: 

dv _ tv 2 — 1 



rfi 1 - t 2 i 



(51) 



where, for simplicity, we have set v = —4-12 1//3 w. Equation ( plj ) has been solved using series expansion by Majorana 
, and this leads to a semi-analytic general solution for the Thomas- Fermi equation (for details see El ) . 



VI. CONCLUSIONS AND OUTLOOK 

In this paper we have generalized a result, derived by Majorana 0], for solving the Thomas- Fermi equation, to 
a wide class of (ordinary) differential equations, that of Majorana scale-invariant equations, as defined in Sect. III. 
We have shown that the search for the parametric solution of such equations of order n can be restricted to that for 
the solution of a differential equation of order n — 1 (and to the computation of one integral involving this solution). 
However, the main result of this paper is not this proposition, which is a direct consequence of known results, but, 
rather, the method and the transformations used to obtain the lower order equation. This has been outlined in Sect. 
Ill and the equation considered has been formally written in Eq. (^3|). In some case this differential equation could 
be highly non-linear, so that in Sect. IV we have developed a variant of the method mentioned above, employing a 
further transformation for the dependent variable involved. Some other simplifications, depending on the particular 
problem considered, can be achieved with a suitable choice for the arbitrary function z{t) present in the parametric 
solution for the differential equation. 

As an illustration, both methods have been applied to reduce the order of Emden-Fowler equations in Sect. V and, 
as a particular case, Thomas-Fermi equation has been considered as well. Remarkably, by using the method of Sect. 
Ill, we have shown that all second-order Emden-Fowler equations can be transformed into first-order Abel equations 
of the first kind. Instead, by using the method of the auxiliary function reported in Sect. IV, the Thomas-Fermi 
equation can be transformed into a suitable first-order equation which can be solved by series expansion. 
We believe that the transformation methods presented here deserve further attention in view of their potential appli- 
cations to scale-invariant differential equations which are of interest for mathematical physics. 
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